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Abstract 

We consider quasilinear, multi- variable, constant coefficient, lattice equations 
defined on the edges of the elementary square of the lattice, modeled after the lat- 
tice modified Boussinesq (ImBSQ) equation, e.g., yz = x — x. These equations 
are classified into three canonical forms and the consequences of their multidimen- 
sional consistency (Consistency-Around-the-Cube, CAC) are derived. One of the 
consequences is a restriction on form of the equation for the z variable, which in 
turn implies further consistency conditions, that are solved. As result we obtain a 
number of integrable multi-component lattice equations, some generalizing ImBSQ. 



1 Introduction 



Recent progress in the theory of integrable difference equation has been quite rapid, 
partly due to the effectiveness of the Consistency-Around-the-Cube (CAC) property as 
an integrability criterion [H El |3]. The work on CAC has so far been mostly restricted 
to one-component equations, but some results do exist for the three-component lattice 
Boussinesq (IBSQ) class of equations [H El El [71 |8]. By this class we mean equations 
relating variables defined on the corners of the elementary square of the cartesian lattice, 
with some equations residing on the edges and some on the square. (In the following 
we will denote the dependent variables by x, y, z and the independent variables n, m, fc, 
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furthermore shifts in the m-variable will be denoted by a tilde, in the n-variable by a hat 
and in the fc- variable by a bar, for example x = Xn+i,m,k+i-) 
One example is the IBSQ (see equations (9) in [7J) 

y = XX — z, y = xx — z, (1.1a) 
z = XX — y + — — —. (l-lb) 

X — X 

Equations f ll.lap are defined on the sides, and the set can be completed by suitable shifts, 
while fll.lbp is defined on the square. Similarly the lattice modified BSQ/Schwarzian BSQ 
(see equations (5.3.19-20) in ^ and (4.7-8) in |5]) is defined by 



yz = x — x, yz = x — x, (1.2a) 
.{zyp^ -zyq^). (1.2b) 



z — z 

(The second equations ( ll.lblll.2bj) can be written in other forms by replacing some shifted 
variables with others, using the first equations.) 

If one eliminates two of the variables and expresses the equation in terms of the third 
one, the final equation turns out to be defined on a 3 x 3 stencil on the discrete 2D- 
plane. In fact, the lattice modified and Schwarzian BSQ equations differ only in their 
one-component 3x3 formulation: ImBSQ ((4.9) in [5j) is obtained after eliminating x,z, 
and ISBSQ ((4.1) in [5]) after eliminating y,zu Further examples of this type are given 
in [8]. 

From the point of view of the CAC approach we can take the quasilinear equations as 
being given on all 12 edges of the consistency cube. In the case of (ll.2ap they are 



yz = 


X — 


X, 


yz 


= X — 


X, 


yz = 


X — 


X, 


(1.3a) 


yz = 


X — 




yz 


= X — 




yz = 


X — 




(1.3b) 


yz = 


X — 


X, 


yz 


= x — 




yz = 


X — 


X, 


(1.3c) 


yz = 


X — 


X, 


yz 


= x — 




yz = 


X — 


X. 


(1.3d) 



Equations fll.3ap should be considered as constraints on the initial data: x, y, y, y, y, z, z, z 
may be arbitrary, but x,x,x are restricted by fll.3ap . Equations ( ll.3b|1.3cp can next be 
used to determine the doubly shifted x, y quantities, resulting with the face equations 

x= ^ ^ , y = -^^^, 1.4 

z — z z — z 

and the corresponding formulae for x, x, y, y. Then we can solve for yx from (ll.3dp and 
the 3D consistency implies, through the third equation, a constraint on z: 

I(z-z) ^ z(z-z) ^ z(z-z) _^ ^^^^ 
z^ 'zz 



*From this point of view the edge equations are considered as Cole-Hopf-hke transformations between 
the main equations, ImBSQ and ISBSQ. 
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This can be separated and we find that 



z z 



(1.6) 



z = 



z — z 



where the F's are some free functions of the parameters and shifted variables of the 
direction indicated by the subscript, i.e., 

Fp = (f){x,y,z,x,y,z,p), Fg = (f){x,y, z,x,y,z, q), Fr = (f){x,y, z,x,y,z,r). (1.7) 

(Here some of the shifted variables can be eliminated by using ( ll.3ap .) 

Up to this point we have only used the quasilinear equations and their consequences. 
The 3D consistency of ;z-shifts, as derived from f 1 1.6 1) and the corresponding other shifts, 
implies further requirements on the (p, which in the case of ImBSQ/lSBSQ are solved by 



In this paper we formalize and generalize the above point of view of starting with quasi- 
linear equations on the sides of the consistency cube and then deriving multi-dimensionally 
consistent multi-component equations. We start in Section 2 by classifying the quasilin- 
ear part (under the constant coefficient assumption) and find that there are only three 
essentially different cases to consider. In Section 3 we derive for each case the consistency 
equations for z and finally in Section 4 solve these equations. 

2 The quasilinear part 
2.1 Classification 

The following properties will be required of the quasilinear equations: 

1. The equations are linear in both the shifted and the unshifted variables. 

2. Precisely two shifted variables appear, we call them x and y. 

3. The third (non-shifted) variable z appears in a nontrivial way. 

4. The coefficients of the quasilinear equations are all constants, or can be transformed 
to constants by suitable [n, m, fc-dependent) transformations. 

The second property is modeled after the existing lattice BSQ results, while the fourth 
property is a technical restriction. It is possible that other integrable models can be found 
if either or both of these conditions are relaxed. 

Clearly the most general equation satisfying the above requirements has the form 



(c.f., (O) 



cj)(x,y,z,x,y,z,p) = —p^. 



(1.8) 
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In order to find the canonical forms we liave at our disposal rotations in x, afiine 
transformations in each coordinate, and z 0\X + + o-^z + t^. Furthermore we can 
apply rational linear transformations of the type 

(X, Z) I— 7- ^, 

where P, Q, i?, S* are linear functions in x, Finally we may sometimes use n, m, /c de- 
pendent transformations, provided that the final result does not explicitly depend on 
n, m, k. 

In Appendix A it is shown that all equations of the above type can be transformed 
into one of the following three forms: 

xz = y \ (A) 

X X = y + z, (B) 

yz = X — X. (C) 

The convention here is that the x- variable appears both with and without a shift, y only 
with a shift and z only without. There is a fourth case of this type, namely x + xy = z, 
but it can be transformed to (A) using {x,y,z) h-> {l/x,y/x,z/x). 

2.2 Reversal symmetry. 

The starting quasilinear equations have a reversal symmetry. Take for example Case B, 
we have 

~ ~ rev shift ^ _ z<r^y ^ ^ 

XX = y + z — y XX = y + z — > xx = y + z — > xx = y + z. 

Thus the reversal must be accompanied with the exchange z ^ y. Here we cannot 
yet determine whether x or the spectral parameter should also change sign, but for the 
solutions found later on it turns out that x should change sign but spectral parameters 
should not. 

Similarly for Case C we have 

^ ^ rev shift ^ ^ y<r^—z ^ ^ 

yz = X — X — y yz = X — X — y yz = x — x — y yz = x — x, 

Here, instead of {x,y, z) (-)■ (x, —z,y) we could have other sign choices. 
Case A is somewhat different, for first we have 

xz = y ^ X — J- xz = y + X — y xz = y + x, 
after which we must make the slightly more complicated transformation to get (jX]), namely 

(x,?/, z) H> (1/x, -z/x, -y/x). 
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Figure 1: Possible definitions of initial values and evolution. Here W stands for the vector 
{x,y,z). When equations are given on the links connecting lattice points some of them 
imply constraints on the initial values while some give part of the evolution. 



2.3 Evolution 

For an equation defined on the elementary square of the cartesian lattice one usually 
considers initial values given on a staircase-like configuration, see Figure [TJ We can, for 
example, have the initial values given on a staircase extending from lower right to upper 
left and with evolution towards upper right, as in Figure [1^). On the elementary square it 
means that the values with no shift or one shift are given as initial values and the double 
shifted quantities are to be derived from the equations. 

In the present quasilinear case, however, we have four equations defined on the edges 
of the elementary square and two of the equations only involve the initial values; these 
equations are therefore to be considered as constraints on the initial data. Thus we can 
take, for example, x, y, z, y, z, y, z as given and then determine x using two of the four 
edge equations, after which the remaining two edge equations give x, y^ for example as 
in f ll.4p . If, on the other hand, we were to consider an initial staircase extending from 
lower left to upper right and evolution towards lower right, as in Figure [T]d), we could 
take X, y^ z, y, z, y, z as given, x, x constrained by two edge equations, and then get from 
the remaining two edge equations two relations involving part of the evolution 

towards lower right. However, in the present case the edge equations are not symmetric, 
since they involve y and 2; in a different way, and therefore the equations giving evolution 
towards lower right and those evolving towards upper right have different structure. In 
the following we will only consider the evolution of Figure [T^) and its completion to the 
consistency cube. 

So far the missing equation in the evolution is the one for z, its form will be determined 
in the following from considerations of multidimensional consistency. 
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3 CAC 



We will now derive the consequences of multidimensional consistency for the basis equa- 
tions. 



Case A 

Writing out the starting quasilinear equations (jSl) we have 

xz = y + x, xz = y + x, xz = y + x. (3.1) 

Now taking various shifts of these equations and then solving for the doubly shifted 
quantities we get 

x=- — y=— — (3.2) 

z — z z — z 

and similar equations for the bar-tilde and bar-hat shifts. There are many other ways 
to write these equations, by inserting suitable combinations of the edge equation. For 
example, a reversal symmetric form is obtained with 

1~_ {x-xy z^_ {xz-xz){y -y) 

X (z — z) {xy — xy) ' x (z — {xy — xy) 

The requirement of 3D consistency on x is derived from 

[xy= (xy= (^)~ 

This and the corresponding equations for y lead to one single condition 

f(i'- ?) +f(^- z) +i(^- F) = 0. (3.4) 

Each of the three terms involve two of the three possible shifts and we can therefore 
separate the variables with the solution 

t=^(Fp-F,), (3.5) 

and similarly (cyclically) for the other shifts. Here Fp,Fq, F^ are some (rational) functions 
associated with the 171,72, k directions (i.e., tilde, hat and bar shifts), respectively, as 
described in (11. 7p . In practice we will eliminate the possible y,y,y dependency using 
(13. ip . The dependency of on the shifted quantities must be nontrivial so that for 
example z ^z. In particular Fp = f3{x, y, z)z is not acceptable. 

Finally the requirement of multidimensional consistency of (13. 5 p implies that the func- 
tion F just introduced must satisfy 

F — F F— F F — F 

'z — 'z ^ —z 1 —z ^ 



t Actually the requirement factorizes, with the other factor being x{z ~ z) + x{z — z) + x{z — z), but 
since this is a first order condition we will not use it. The same happens in the other cases. 



6 



where Fp = ^, x, z, p), etc. This is the main equation to be solved in the next 

chapter. Note again that the spurious solution Fp = (3{x, y, z)z is not acceptable. In 
solving this equation we should use the previous equations f l3.2p .f l33|) to eliminate all 
doubly shifted quantities, but we also have to eliminate shifted y's using (13. ip . 



Case B 

The equations are 



xx = y + z, xx = y + z, xx = y + z. (3.7) 



After computing the various shifts of these equations we use them again to eliminate 
shifted y's from the results. The doubly shifted variables then have the form 

c: z — z' ^ xz — xz: 

X = y = ^ ^ . 3.8 

X — X X — X 

These equations are not symmetric under the reversal, but it is easy to add linear combi- 
nations of (13. 7p to get reversal symmetric forms, for example of the form 

{x — x){x — x) = {z — y) — {'z — y), {y — z){x — x) = (z — y)x — {'z ~ y)x. (3.9) 

In either case the consistency condition becomes 

z{x ~x) +^{x — x) + z{x — x) = 0, (3.10) 

which is resolved by 

t=:^(Fp-F,). (3.11) 
X — X 

The 3D-consistency for z then leads to 



F — F F — F F — F 



(3.12) 



This is similar to the result (13. 6 p for Case A, but note the different denominators, x — x 
vs. z — 

Case C 

The starting equations are now given by 

yz = x — x, yz = x — x, yz = x — x. (3.13) 
Following the procedure given above we first get 

^ xz — xz' ^ x — x y — y 

X = — y = —— — - = —z— — -. (3.14) 

z — z z — z z — z 
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Again it is easy to generate the reversal symmetric forms, for example, 



x{z — 2)+ x{y — y) = {y + 'z)x — {y + ^)x, y{z — z) + z{y — y) = 2{x — x). (3.15) 

or _ _ ^ ^ ^ 

{x-x){y-y^ + z-z) ^ 

{y-y){z-z) ' {y-y){z-z) 
From any form we get the consistency condition 

z(z — z)z + 2(2 — z)z + z(z — z)2 = 0. (3-17) 

The form of z must then be 



^ zz 

z = .^-^{Fp - Fg) (3.18) 



and its 3D consistency implies 



F — F F — F F — F 
F - F ~ F - F ~ F - F ' 



This is essentially different from Case A (13. 6 p or Case B (I3.12p . 



(3.19) 



4 Solutions 

We will now proceed to solve the final CAC equations. They are all of the form 



F - F F - F F - F 



p \J \J q \J p \J y \J q 



(4.1) 



where the subscript indicates on which variables the function depends on, c.f., (II. 7p . 
The meaning of G depends on the particular case. In the shifted functions, for example 
Fp = (j){x, y, S, X, z, p), we should use the previous formulas to express x and z in terms of 
the once-shifted quantities, and thus we should think Fp as a function of 

Equation (14.11) contains functions of at most two of the three sets of shifted variables 
appearing in the equation. Concentrating on p, q sets, say, (and fixing the values of 
r, X, y, z) we find easily that Fp must have the form 

Fp = {Gp - C) + 5,, (4.2) 

Kjp — Kjq 

where A, B are some functions of the unshifted variables, and of the shifted variables 
and parameters indicated by the subscript, while C may only depend on the unshifted 
variables. Using the same form for the other F's as well solves the equations (14. ip . Thus 
the final problem is to solve (14. 2p . 
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Case A 



Using ffOia. 203.51) the explicit form of is 

,/ x-x Fp-Fg \ A{x,z,p) - A{x,z,q) 

[x,y,z,- — — —,p = [z-C) — + B{x,z,q) (4.3) 

\ z — z z — z J z — z 

where the dependence on unshifted variables is hidden. Let us first consider the generic 
case in which y4(5;, i) is more complicated than cz. Then as a function of x,z the RHS 
has a simple pole at ? = z . Observing the form of the LHS we see that if Fp — Fq gtx — x 
(i.e., z x) (f) must be linear in the fourth and fifth variables, otherwise it would also 
have other poles. Thus we find that 

Fp = apX + f3pZ + 7p, (4.4) 

where the functions ap, f3p, 'jp may depend on x, y, z,p (but not on shifted quantities), and 
similarly for Fq,Fr. This even includes the special cases A = cz and Fp — Fg ocx — x. 
Equation f l4.3p can now be written in the form 

ap{x -x) + Pp{apX+/3pZ + 7p - a^x - (3gZ - -fq) + (z - z)% 

= {z- C){A{x,z,p) ~ A{x,z,q)) + B{x,z,q){z -z). (4.5) 

Concentrating on the tilde-shifted quantities we find that 

~ N ci{p)x + C2{p)z + csip) 

A{x,z,p) = — , 4.6 

z — C 

where the functions Ci{p) may depend on x,y,z, and similarly for the hat-shifted quanti- 
ties. When these are used in (14. 5 p the coefficient of a; factorizes as {'z—C)[ap+apf3p—Ci{p)]. 
In explicit form the equation to solve is 

a{x, xz — X, 'z,p) + a{x, y, z,p)f3{x, xz — x, z^p) = ci{x, y, z,p), (4.7) 

where we have already eliminated fusing (13. ip . Equation (14. 7p is solved in Appendix lB.il 
Once the solution for a, /3 is known it is relatively straightforward to solve the remaining 
equations. 

There are three different cases: 



Case A.l: For a = the remaining equations imply /5 = 6o + -2, 7 = 7(p), leading to 
the solution 

^=bo + z + (A-1) 
z — z 

where we have redefined the spectral parameters 7(p) — i- p etc. Here bo can be eliminated 
by the transformation Znm ^ ^nm + |(?t- + m + A;)6o5 but in order to keep intact the form 
(13. ip this needs to be accompanied with y ^ y \ hQ\{n \ m \ k)x. 

Equation (lA-ip depends only on z so the system is triangular with z driving the 
system, x determined next from (13. 2p and y last from (13. ip . The z-equation is that of the 
discrete potential KdV: (z — z){z — z") = p — q. 
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Case A. 2: For a = a{p) 7^ 0, /3 = (3{p) the remaining equations imply that a, (3 are in 
fact constants and 7 does not depend on lattice parameter. This leads to a trivial solution 

z = bo + aox, i.e., z = bo + a^x, 

which in turn implies that x is constant. 



Case A. 3: For the third case, starting with a{x,y, z,p) = a{p)/x, P{x,y,z,p) = y/x + 
b{p), the conditions from the remaining equations imply that b{p) is a constant and that 
7 does not depend on lattice parameter. This leads to a genuine coupled system with z 
given by 

^ , , I px — qx /.X 
z-y x = bo + - ^ I ■ (A-2) 
X z — z 

The parameter 60 can be eliminated with Znm ^ Znm + bo\{n + m + k) , y ^ y + bo\{n + 
m + k)x, and we have redefined a{p) — )• p. This equation is symmetric under the reversal 
if we replace the divisor x in the last term using x = (xy — xy) /(x — x): 

~ y/^ _ {px - qx){x - x) 
(z -z)(xy-xy) ' 

We recall that in this case the accompanying equations are (13. ip together with (13. 2p or 
(13. 3p . Furthermore, the triply shifted quantities are 

i = X f/f'^.^+'if'^+'g'^l, , (4.8a) 
px[z — z) + qx{z — z) + rx[z — z) 

S US px(x — x) + qxix — x)+ rxix — x) , . 

X px[z — z) + qx{z — z) + rx[z — z) 

f = z-x i;L^'^^+i^f"^+'^!"^l, . (4.8c) 

px[z — z) + qx{z — z) + rx[z — z) 

Case B 

This case is similar to Case A, differing in the exchange of x <H- ^ in some parts of the 
formulae. For example the equation corresponding to (14.30 now reads 

,/ z-z Fp-Fq \ A{x,z,p) - A{x,z,q) j,^^ ^ . 

(p{x,y,z,- — — ^,p = [x-C) — + B{x,z,q), (4.9) 

\ X — X X — X J X — X 

and as before we conclude that Fp must have the form (14.41) . When this is used in (14.90 
we get, corresponding to (14.50 

ap{z -z) + Pp{apX+(5pZ + 7p - UqX - PgZ - 7^) + {x - x)% 

= {x - C){A{x,z,p) - A{x,z,q)) + B{x,z,q){x -x), (4.10) 
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and then _ _ 

X — O 

The coefficient of z now yields 

a(x, XX — 2,p) + f3{x, y, z,p)f3{x, xx — z, z^p) = C2{x, y, z,p). (4-12) 

This differs from (]4.12p mainly by being quadratic in /3. It is solved in Appendix IB.2I with 
the result f lRT2l) . 

When the solution ( IB. 121) is substituted into fl4.10p the subsequent equations quickly 
imply that the parameters Oq, b^, bi are in fact constants. There are two solutions: 

Case B.l: The ffist solution is obtained if 7 = 0, 61 7^ 0: 

z = aQ — bobi X — b\y + ipQ + bix)x (4-13) 

Recall that the other equations that go with the above are (13. 7p and (13. 8 p or ( 13. 9p . 

Equation f l4.13p can be simplified by further transformations which nevertheless must 
keep the equations fl3.8p intact. We find that there are two canonical forms, the ffist of 
which is 

z = XX — y. (B-1.1) 
If 6^ 7^ 1 we can use the linear transformation 

{■^,y,^)^y-^ b^^i^y -^bi-i ^ (fei+i)(6i-i)2 > ^ bi~i ^ {bi+i){bi~i)^ ^ ^ 

and follow it with the scahng {x,y,z) (s2(n+m,+A:)+i^^ _g4(n+m+fc)^^ ^4{n+m+fc)+4^-j^ where 
= bi, to transform fl4.13p into the form (]B-l.ip . If 61 = —1 we use instead 

(x, y, z) {-x+\bQ, y-x\bQ+{\aQ+\bl){l-n-m-k), z-x\bo+{\aQ+\bl){n+m+k) + \bl), 

and scaling. If 61 = 1, 60 = we get (IB-l.ip directly. 

If 61 = 1, 60 7^ we cannot eliminate 60, only scale it, but we can change the value of 
Oq with n, m, k dependent translation. This yields the second form 

t= (x + 1) (§- 1) -y. (B-1.2) 

Both of the forms ( IB-l.ip and ( 1B-1.2P are manifestly reversal symmetric, but the 
second case requires that x changes sign under reflection. 

The triply shifted quantities reveal that the system (IB-l.ip is in fact periodic on the 
diagonal of the cube: 

X = X, y = y, z = z. (4.14) 



In the other case (1B-1.2P we have a similar result 

X = X + 1, y = y + x + l, z = z + x. (4-15) 
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Case B.2: The other solution is obtained for 7 = p, 61 = 1: 

z = Qq — bo(x — x) — y + XX + — — —, (B-2) 

X — X 

together with fl3.8p . Here the parameter can be transformed away by an n,m,k de- 
pendent translation, similar to the ones used before. For 60 = we recover IBSQ (11.11) . 
This result is reversal symmetric, but if x changes sign, as in the previous results, then 
p, q should not change sign. 

The triply shifted quantities are 

S {q-r)x + {r-p)x + {p-q)x i« n 

X = bo + x + —^ — - , — , — ^, (4.16a) 
x{z — Z) + X[Z — Z) + x{z — Z) 

5 (g-r)?+ {r-p)z + {p - q)z (a Tau\ 

y = box + y+— — , — ^ — ^, (4.16b) 

X[Z — Z) + X[Z — Z) + X[Z — Z) 

~ T S (q' — r)xx + (r — p)xx + (p — q)xx 

z = z + box — — — — — (4.16c) 

x{z — z) + x[z — z) + X[Z — z) 

Case C 

For this case we have G{p) = Fp in (14. 2p . i.e., 

y-yzz{Fp-Fq) \ A{y,z,p) - A{y,z,q) 

x,y,z,-z- — -, — ,p ={Fp-C) + B{y,z,q), 4.17 

z — z z — z J Fp — Fg 

It is now more convenient to eliminate shifted x-variables, since they appear linearly, and 
therefore we use (l){x,y, z,y,'z). In the generic case the pole structure now implies 

Fp = [apy + /3pl + jp]/I, (4.18) 

and then (I4.17P becomes 

z{y-y)ap + ^{apyz-aqyz + (/3p - (3q)'zz + -fpZ - -fqz) + {z-z)yp 
= z{apy + l3pZ + -fp - Cz) [A{y, z, p) - A{y, z, q)] + 

B{y,z,q)[apyz-aqyz + {(5p - (5q)zz + ^pZ - ^qz\. (4.19) 

Since the dependence on tilde-shifted quantities is still linear, except for the A term, we 
conclude _ _ 

A(y.J.p) = £lM|l:|<flI±^. (420) 

Upy + [3pZ + 7p - Cz 

When this is substituted into (14.170 the result is linear in y, 'z and a suitable linear 
combination of their coefficients yields the equation 

a{yz + x,y,z,p) (z/z) [a{x,y, z,p)y + -f{x,y, z,p)] - -/{yz + x,y,z,p) a{x,y,z,p) = 

a{x, y, z, p)cs{x, y, z, p) - 7(0;, y, z, p)ci{x, y, z, p). (4.21) 

This equation is solved in Appendix IB. 3 1 with four solutions, a = 0, (IB.lSp . (IB.28p . and 
(IB.Sip . The other equations for each case are (I3.13P and (I3.14p or (I3.15p . 



12 



Case C.l: a = 0. First assume that also 7 = 0, then the remaining equations can be 
easily solved with the result 

t=^3^. (c-i.i) 

z — z 

Next if 7 ^ the analysis can be divided into sub-cases according to whether or not 
dzl = 0. In the former case only trivial result follow [z oc z), in the latter case one gets 
the solution 

z — z 

Both of these results depend only on z and its shifts so they form a triangular system, as 
in Case A-1. In the case of (1C-1.2P the equation can be scaled to lattice mKdV. 

Case C.2: 'y{x,y,z,p) = a{x,y, z,p){x + g{p))/ z. Concentrating next on the part 
multiplying y and absorbing all q dependency into B there remains 

-13 (x, y, z, p) + a{x, y, z, p)/a{x, y, z, p) + a (x, y, z, p)/a{x, y, z, p) + B'{x, y, z, x, y, z) = 0. 

This equation is similar to f l4.7p but has the additional B' term. Thus the solution is also 
a bit more general, 

a{x,y,z,p) = ai{p)a{x,y,z), l3{x,y,z,p) = l3i{p) + b{x,y,z), 

with suitable values for the auxiliary functions ci,B'. Form the coefficient of z we first 
find that is a constant and then get the equation 

a'ix,y,z)/z'^[ai{p){g{p) + x) - ai{q){g{q) + x)] = C2ix,y,z,p) - C2(x, z, g). 

There are now various way to solve this. 

1: If a{x,y,z) = z'^a^ then ai{p) must be a constant, which can be absorbed into Oq. 
We can take g{p) cc p and then we get a solution with an extra term in comparison with 

t=_^~+;,^iZ^. (C-2.1) 

z — z 

The triply shifted quantities are 

I ^ Xjqz - rz) + X{rz - pz) + X{pz - qz) ('4 22a) 

^ X{z-^)+X{z-z) + X{z-z) ' 

S ^ 1 ^ 1 J(9 -r) + z{r- - p) + z{p - q) (4 22b) 

z zx(z-z)+X{z-z)+Xiz-zy 



a Xzz{q — r) + Xzz{r — p) + Xz'zijp — q) 
X{z-^)+X{z-z) + X{z-z) ' 

where we have used the shorthand notation X = x — p, X = x — q, X = x — r. 



(4.22c) 
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2: If a{x, y, z) = z'^a^/x then we must have g{jp) = aQ/ai{p). After choosing ai{p) oc p 
and scahng this yields the solution 

t=-cio- + -^^t^. (C-2.2) 

X X z — z 

If (io = we can divide this equation with the x equation and write the result in terms of 
w := zjx in the form of lattice KdV. 
For triply shifted quantities we get 

S , xxzin — o) + xxJ(q — r) + xx'zir — p) , . ^ 

X = -do-y ' , ' 4.23a 

X px{z — z) + qx[z — z) + rx[z — z) 



X x[z — 'z)+ x[z — z) + x[z — z) 
z px{z — z) + qx(z — z) + rx{z — T) ' 
pJiz^ —z) + qz(z — z)+ rz{z — z) 



(4.23b) 



px{z — 'z) + qx{z — z) + rx{z — z') 
pxzz{q — r) + qxzz{r — p) + rxz'z^p — q) 
px^z — z) + qxiz — z) + rx{z — 'z) 

Both ( lC-2.11) and (lC-2.21) provide a closed system in x, z so they are really two com- 
ponent systems. 

3: If a{x,y, z) does not have either of the other forms then both ai{p) and g{p) must 
be constants. This leads to a trivial solution: One finds first that zz = —a{x, y, z){x + gi) 
but then the third shift yields x + gi = 0, z = 0. 

Case C.3: a{x,y,z,p) = a{p)z/y, 7(x,?/,2;,p) = a{p)[gi{p)x + g2ip)]/y, gi ^ 1. Since 
the X, y, 2;-dependency is now explicit it is rather easy to solve the remaining equations. By 
considering different powers of y,^ one finds that gi{p) = di/a{p), that /3 cannot depend 
on the spectral parameter and can therefore be eliminated, and that g2{p) = d2/a{p). One 
also finds values for the auxiliary quantities ci, C2, C3, B. The resulting equation is 

- d2X + di zpyz-qyz 

z = \ — - — . (C-3) 

y y z- z 

If (ij = we recover the mBSQ/SBSQ equation (11.21) . In this case we get a bit simpler 
results for the triply-shifted quantities if we eliminate the shifted x-variables. 

S _ d2X + di^ zyy{q - r) + zyy{r - p) + zyy{p - q) 

y z[qy - ry) + z[ry - py) + z[py - qy) 

I = y ^ z{y-V) + z{y_-y)+ny-y) .4_24b) 

z{qy - ry) + z{ry - py) + z{py - qy) ' 

~ , /, 1 \ (q - r)z + (r - p)z + (p - q)z n 
z = d2Z + {di + d2x)—^ — zr, — 4 - ~n ~ — ^ (4.24c) 

z{qy - ry) + z[ry - py) + z[py - qy) 

^ ^^ qrzjy -y)+ rpzjy -y)+ pqz{y - y) ^ 

z{qy - ry) + z{ry ~ py) + z{py - 
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Case C.4: a{x,y,z,p) = a{p)[x + g3{p)]z/y, -f{x,y,z,p) = a{p)[g2{p) + {x + g3{p)){x + 
gi{p))]/y, g2 7^ 0. This is similar to Case C.3. One first finds that (3 = 0, then that a{p) is 
a constant, which can be scaled to —1, (73 (p) = do — gi{p) and (72 (p) = QiipY — dodiip) + d'l- 
Since p-dependency enters only through gi{p) we put gi{p) = p + do/2, furthermore by 
^-translation we can eliminate do- The result is therefore 

- xx + di zpyz-qyz 

z = \ — (C-4) 

y y z- z 

This can be written in a reversal symmetric form, if we replace the z on the RHS using 
z = (x — x)/ {y — y), which yields 

% = x2 + d,+ (C-4r) 

{y-y){z-z) 

Furthermore: 

X = diy/y+[z{q{x — x)x — r{x — x)x) + z'{r{x — x)x — p{x — x)x) + 

z{p{x — x)x — q(x — x)x)]/ A, (4.25a) 

y = y[^{x — x) + z'{x — x) + z(x — x)]/A, (4.25b) 

z = {diz['z{x + q — X — r) + ^(x + r — x — p) + z{x + p — x — q)] 
+z[z{qx(x + r) — rx(x + q)) + z{r-x{x + p) — px{x + r)) + 

'z{px{x + q) — qx{x + p))]} I A, (4.25c) 
where A = i'[(g — x) (x — x) — (r — x) {x — x)] + 
z[{r ~ x){x — x) — {p — x){x — x)] + 

z[{p — x)(x — x) — (g — x)(x — x)]. (4.25d) 



5 Discussion 

We have considered multi-dimensionally consistent equations generated by quasilinear 
equations of three variables defined on the edges of the consistency cube. Assuming 
constant coefficients and dependency on only two of the three possible shifted variables 
these equations were first reduced into three canonical cases (jX]), ([B]), or ([U]). In each 
case, starting with the equations given on the 12 edges of the consistency cube, we can 
propagate x to other 7 corners of the cube, which leaves 5 equations. Assuming free 

initial values for y, y, y, y we can use 4 of the remaining equations to obtain y, y, y, y 
and the final equation gives the first consistency requirement for z, ( 13. 4p . f IS.lOp . or fl3.17p . 
respectively. The equations for the shifted z's are equations on the sides of the cube with 
additional consistency conditions, this implies equations for the remaining freedom in the 
F function, given in (13. 6p . (I3.12p . (I3.19p . respectively. These functional equations were 
solved with the results given above. The results can be summarized as follows: 
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One component results. In three cases the resulting z equation only involved z- 
variables, see ( ]A-ip (KdV), ( 1C-1.1|) . and f lC-1.2p (mKdV). For these equations one can 



reverse the point of view and state that for these z-equations one can adjoin edge equations 
(of the canonical form) consistently. One may ask whether any consistent 1-component 
equation (such as those in the ABS list [3]) can be so extended, if more freedom is given 
for the form of the edge equations. 

Quasilinear z-equation. In the case of f l4.13p . or its two canonical forms ( IB-l.ll) and 
( IB-1.2p . the z equations are also quasilinear, and in fact reversal symmetric. However, 
these equations do not allow evolution towards lower right or upper left, because then 
these equations do in fact impose further conditions on the initial values. The triply 
shifted quantities f l4.14p or (14.151) suggest periodic reduction along the diagonal. 

Two component results. For Case C we found two essentially 2-component results 

z = —zx + z— — —, (^-2.1) 

z — z 

and _ ^ 

I=-do'- + ^ P^l-r , (C-2.2) 
X X z — z 

since together with the x-equation in (I3.14p the variables x, z form a closed system. How- 
ever, since neither x nor (lC-2.ip involves x, this system does not define evolution in all 
directions. If do = ( 1C-2.2P can be reduced to ImKdV for w := z/x. 

Fully three component results. We found four fully three-component systems that 
allow evolution in all directions. One criterion for this is that the z equation must involve 
un-shifted y, which does not appear in any of the edge equations. 

• For the base-equation (jA]), i.e., xz = y + x, we found the 2;-equation, 

t=^ + l^!-r, (A-2) 

X X z — z 



or in the reversal symmetric form flA-2rP 



• For the base equation ([B]), i.e., xx = y + z, we found a mild generalization of the 
IBSQ equation ffTTTj) : 

z + y = Bq^x — x) + X X + — — (B-2) 

X — X 

• For the C-case, i.e., zy = x — x, we found two equations, which are kind of modifi- 
cations of the ImBSQ/lSBSQ equation (11.20 : 

- d2X + di zpyz-qyz 

z = \ — - — , (C-3) 

y y -' 

and 



pyz 


- qyz 


z 


— z 


pyz 


- qyz 



+ - ^^: . (c-4) 

y z- z 
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A Classification of the quasi-linear part 

Introducing vector notation 



u={x,y,l) , U = {x,y,z,l) , 

we can write the equation as 




(A.1) 



The rotation freedom means that we can introduce primed quantities by 



u = mu , 




Si2 

S22 ty 

02 03 tz 

1/ 



(A.2) 



and this induces the transformation A ^ A' = m^ AM. 
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A.l Classification of the quadratic part 

Let us first consider the quadratic part, which is determined by 



^ an ai2 a,,, 

V«21 022 023/ 



and which can be transformed as Ag ^ A' = m* Aq Mq by 



Sll Si2 
S21 S22 
Oi 02 03, 



= ( III '22 J ' = I ° I • ^^-'^ 



The main classification is as follows: 
• If one or both of ai3, 023 are nonzero then we can choose Sij so that aissn + 023521 = while 
013S12 + 023S22 7^ and S11S22 — S12S21 7^ 0. Then using 01,02,03 we can put A'^ into the form 



f 9 
1 



This we divide into three cases: 

Case 1: f = g = after which we exchange x, y and get 

xz = (linear terms). 

Case 2: / = 0,5 7^ 

xy + yz = (linear terms). 
Case 3: Finally if / 7^ we can use the freedom still left in Sij to put 5 = 0: 

XX + yz = (linear terms) . 

• Next if ai3 = 023 = then oi , 02 are superfluous and we only have the s-part to operate with. 
Let us denote the left 2x2 block of Aq by Ag^ then the remaining transformation is s^AgS. This 
can be used to bring Ag into one of the following 6 normal forms 

0\ /O 1\ (I 0\ (I 0\ (I a\ f 1 

0) \o 0) \o oj'Vo ij'Vo v'V-1 

where the parameter a is nonzero. The first case must be omitted, as it leads to linear equation. 
Thus we have 5 more cases: 
Case 4: 

xy = (linear terms). 

Case 5: 

XX = (linear terms). 

Case 6: 

XX + yy = (linear terms). 

Case 7: 

XX + axy + yy = (linear terms) . 

Case 8: 

xy — yx = (linear terms). 
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A. 2 Detailed classification 

We will next do a detailed analysis and classification of the quasilinear equations. 



Case 1 

In Case 1 we must have 62 7^ in order to have y dependency, therefore we scale to 62 = — 1- 
This case is invariant under the full transformation (IA.4p . and if we take sn = 1, S12 = 0, S22 = 
1, 01=02 = 0, 03 = 1, tz = S21 — bi, tx = — C3 we can simplify the linear terms into 

xz-y+{Bt + A)x + By + s{B - 1) + C = 0. (A.5) 

This is further divided as follows: 

If B = y disappears and we can eliminate the constant using s. Now we must have A 
otherwise the system leads to triviality, thus we can scale to 

xz = y + X, (A. 6) 

which is Case A 

If i? ^ we can simplify the system by eliminating x using t and obtain 

xz = y-By-s{B-l)-C. (A.7) 

li B ^ 1 we can next eliminate the constant C with s and after that redefine (x, y) 1— )• 

(^Qu+m+ky^ Qu+m+k^^ obtaining 

yz = x — X, (A. 8) 

which is Case C. If instead = 1 we can eliminate C with (x, y) ^ {y + {n + m + k)C,x), which 
again yields ([C]l . 

Case 2 

The quadratic part of Case 2 is invariant under (IA.4j) if su = 1, •S21 = 0, S22 = 1- Choosing 
furthermore S12 = — 02, oi = 0, 03 = 1, tz = 621 ty = —bi, tx = 02(61 + ci — C3) — C2 the equation 
becomes 

xy + yz = A'x + B'z + 02(5' - A')A' + C'. 
li B' ^ the transformation 

'xt + y + s B' -A'/B'{xt + y + s)-D'{x + l) + z + t" 



{x,y,z) ^ 



x + l'x + l' X + 1 



takes this to i^AM with B = A'/B', A = C = D' where D' = 02A'{B' - A')/B' + C'/B'. 
If B' = 0, A' ^ then we can eliminate the constant term with 02 and scale to 

xy + yz = x. 

This can be transformed to by {x,y,z) 1— t- (l/y,—x/y,—z/y). 
If A' = B' = then C" 7^ and can be scaled to 1 and we have 

X y + y z = 1. 

With {x,y,z) 1-^ {y/x, l/x,z/x) this can be transformed to the third canonical form 

XX = y + z, 

which is Case B. 
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Case 3 

For Case 3 and the other remaining cases the possible rotational freedom in mq is not useful. 
Using the translational freedom with 03 = l,tz = —b2,tx = —bi,ty = — C3 we can bring the 
equation into the form 

XX + y z = Ax + By + C, (A.IO) 



If C = but AB 7^ then the equation is taken to ([A]) with the transformation 

(x, y, z) ^ {-A{x + l)/{y - 1), {A^/B)x/{y - 1), -B{x + z + l)/(y - 1). 

li C = B = the system is trivializes with (x, y, z) 1-^ (x + A, y, z{x + A)). 
If C = yl = we scale to = 1 and get 

X X + y z = y. (A. 11) 

This is transformed to (jB]) with {x,y,z) 1-^ {x/y, l/y, —z/y): 
Now we assume C ^ and scale it to 1: 

xx + yz = Ax + By + l, (A.12) 

Here, for convenience, we parametrize A = —F + \/F, then we can transform ([02]) to (jXTt) 
(with B = s = 0, C = -1) using 



{x,y,z) ^ 



-(F^y + l) -X -BF^y + z - B{F -1/F) 



Bx + Fy ' Bx + Fy' Bx + Fy 



Case 4 



In this case we must have 62 7^ 0. We use the freedom in the z transform Oj and the translations 
ti to bring the equation into 

xy = y + {Aoi + B)x + Az + (At^ + C). 

Now if j4 = the variable z disappears completely, which is not allowed, thus using oi,tz and 
scaling we obtain 

xy = y + z, 

and with {x,y,z) i->- {—y/x,l/x,—z/x) this becomes 
Case 5 

Here we must have 62C3O3 ^ 0. We can then choose 01,02, tz,tx and scaling so that we obtain 
Case 6 

The situation is similar to Case 5 and leads first to 

X X + yy = z. (A. 13) 

With (x, y, z) 1-^ ((y + l)/x, iy/x, z/x) this gives the B = —1, C = —1 subcase of (lA.Tp . 
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Case 7 

Now using the translation freedom and su = 0, su = —1, S21 = 1, S22 = a we get 

X X + axy + yy = z, a ^ 0. (A. 14) 

For convenience let us parametrize a = i{n — Then with transformation {x,y,z) 1— > 

{{Ky + l)/x,iy/x,Kz/x) this case becomes (IA.7P with B = — C = —k~^. 

Case 8 

The rotational freedom is not useful but we can effectively use translations and get 

X y — yx = z. (A. 15) 

With (x, y, z) ^ (1/y, —x/y^ z/y) this becomes ([C|) . 

In conclusion, by suitable transformations we can bring all nontrivial models satisfying the 
condition mentioned at the beginning of Sec. 2 to one of the 4 forms ([A|), ([C]), or ([B]) . 

B Solving functional equations 

B.l a + ajS = c 

This is the equation needed in Case A, 

a{x, xz — X, ^iP) + a(x, y, z,p) /3{x, xz — x, 'z,p) = ci(x, y, z,p). (B.l) 

It is obvious that this has the solution a = 0. Furthermore, since 

f{x,xz - x,z,p) = g{x,y,z,p) ^ f = f{p), (B.2) 

we also find that if a ^ then a{x,y,z,p) = a{p) ^ f3{x,y,z,p) = I3{p) and this provides 
another solution. 

Next let us assume that both a and /3 depend non-trivially on some of the other possible 
variables. Computing the y-derivative of (jB.ip shows that since f3 must depend on some of its 
first three variables a cannot depend on its second variable. Thus we get 

a(x, z,p) + a{x, z,p) I3{x, xz — x, z,p) = C2{x, z,p) (B.3) 

Now by operating on this with {xdx + dz)dz we obtain 

[xdxa{x, z,p) + dza{x, z,p)] dzl3{x, xz — x,z,p) = 

If the first factor vanishes then we are back to a = a{p) and therefore we must assume that 
/3{x,y, z,p) = P{x,y,p). But then it follows that a{x,z,p) = a{x,p) and the equation reduces 
to 

a{x,p) + a{x,p) /3(x, xz — x,p) = C2(x, z,p) (B.4) 
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Now computing the z-derivative and then applying ()B.2p we find that f3{x,y,z,p) = {y/x)hi{p) + 
60(2;, p) and then the equation becomes, after some rearrangements, 

a{x,p) + a{x,p)[-{x/x)hi{p) + hQ{x,p)] =C2{x,z,p) (B.5) 

Now it is clear that a(x, y, z^p) = ai{p)/x + ao{p) and bQ{x,p) = h2{p)/x + h^{p). Next we find 
that since ai / we must have 62 = = 0, 61 = 1. 

Thus collecting the results we find that equation (jB.ip has precisely the following solutions: 

• a{x,y,z,p) = 0, 

• a{x,y,z,p) = a{p), (3{x,y,z,p) = P{p), 

• a{x,y,z,p) = a{p)/x, I3{x,y,z,p) = y/x + (5{p). 

B.2 S + /3^=c 

The equation to solve for Case B is 

a(x, XX — z, z,p) + fi{x, y, z,p)/3{x, xx — z, z,p) = C3(x, y, z,p). (B.6) 

Again we have the solutions a = a{p), (3 = /3{p) and next we assume that there is some other 
dependency as well. 

Computing the y-derivative and recalling (IB.2P we find that /3 cannot depend on its second 
variable and the equation becomes 

a(x,xx - z,z,p) + f3{x,z,p)f3{x,z,p) = cz{x,y,z,p). (B.7) 

Now dividing both sides with /3(x, and then operating with {dx + xdz)d^ we get 

= dza{x,xx- z,z,p) [{dx + xdz)l/l3{x,z,p)] . (B.8) 

Since /3 cannot just depend on p the first factor must vanish and thus a cannot depend on z, 
then from (|B.7p we conclude that neither can /3. Thus we find 

a{x,xx - z,p) + /3{x,p)l3{x,p) = C3{x,y,z,p). (B.9) 

Next computing the z-derivative of this we find that a{x,y, z,p) = yai{p) + a2{x,p), thus we 
get 

xxai{p) + a2(x,p) + f3{x,p)f3{x,p) = c!^{x,y,z,p). (B.IO) 

Since j3 has to depend on its first variable this dependency must be linear, /3(x, y, z,p) = xbi{p) + 
60 (p), leading to 

xxai{p) + a2{x,p) + [xbiip) + bo{p)][xbi{p) + bo{p)] = c!^{x,y,z,p). (B.ll) 

Now also a2{x,p) has to be linear and the result is therefore 

{ a{x,y,z,p) = ao{p) - xbo{p)bi{p) - ybi{pf (B 12) 

1 f3{x,y,z,p) = bo{p) + xbi{p) 

With the special case bi{p) = we get the solution mentioned at the beginning of this 
section. 
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B.3 a{z/^){ya + 7) — 70^ = c^^a — C17 

The equation to solve for Case C is 

a{yz + X, y, z, p) (z/z) [a{x, y, z, p)y + 7(3;, y, z,p)] - -f{yz + x, y, z, p) a(x, y,z, p) = 

-q(x, y, z, p)c3{x, y, z,p) + 7(x, y, z,p)ci (x, y, z, p). (B.13) 

Clearly this has the solution a = 0. Next we assume a ^ 0. After operating on ()B.13P divided 
by a(x, y, z,p) by dyd^ the result factorizes as 

[a(x, y, z,p) dy-f{x, y, z,p) - j{x, y, z,p) dya{x, y, z,p)] x 

[a{yz + x,y,z,p) - zd^a{zy + x,y,z,p)] = (B.14) 

a) Let us first consider the first factor and assume that the second factor does not vanish. 
Then clearly 

l{x,y,z,p) = a{x,y,z,p)g{x,z,p). (B.15) 
Now operating on ()B.13P divided by a{x,y,z,p) by {dz — ydx)dz the result factorizes as 

[zydxg{x,z,p) - zdzg{x,z,p) - g{x,z,p) -y]x 

[a{yz + x,y,z,p) - zd^aiyz + x,y,z,p)] = (B.16) 

Again we use the first factor, and since g does not depend on y we get two equations 

zdzg{x, z, p) + g{x, z, p) = 0, zdxgix, z,p) -1 = 0. 

which are solved by 

g{x, z,p) = (x + g{p))/z. (B.17) 

This yields the solution 

a = arbitrary, 'y{x,y, z,p) = a{x,y, z,p){x + g{p))/ z. (B.18) 

b) Next we take the second factor in (jB.14p . which is solved by 

a{x,y,z,p) = za{x,y,p). (B.19) 
Then operating by on (IB.lSp yields 

a{x,y,p)d^g{yz + x,y,z,p) = 0, (B.20) 

so that 

l{x,y,z,p) = g{x,y,p). (B.21) 
The z derivative of (lBT3]l now vanishes and next we consider its y derivative, which yields 

dya{yz + X, y,p)[zya{x, y,p) + g(x, y,p)] 

-dgg{yz + x,y,p)a{x,y,p) + z a{zy + x,y,p)a{x,y,p) = 0. (B.22) 
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Since z-dependency is now explicit it is useful to change variables by replacing z with x using 
= X — X, which also implies dy ^ dy + zd$. Equation (IB.22P then becomes 

[ydya{x, y, p) + {x - x)d^a{tj, x,p)][x - x + g{x,y,p)/a{x,y,p)] 

-[ydy9{x,y,p) + {x - x)dri:g{y,x,p)] + {x - x)a{y,x,p) = 0. (B.23) 

Now computing the y-derivative of this yields 

[ydga{x, y, p) + {x - x)d^a{y, x, p)] [a{x, y, p)dyg{x, y, p) - g{x, y, p)dya{x, y, p)] = 0. (B.24) 

The first factor vanishes only if a{x,y,p) = a{p) and then (IB.23P implies 

ydyg{x, y,p) + x[dsg{y, x,p) - a{p)] = 0, dj,g{x, y,p) - a{p) = 0, (B.25) 

which is solved by 

9{x,y,p) = a{p){x + g{p)). 

This yields, however, only a subcase of (jB.lSp . 

For the second factor in (]B.24p to vanish we take 

9{x,y,p) = a{x,y,p)gi{x,p). (B.26) 

Substituting this into (|B.23p and then operating with 5^ yields 

yd^a{x,y,p)dlg{x,p) - d^a{x,y,p)[{x - x)dlg{x,p) + 2dxg{x,p) - 2] = (B.27) 

We divide solving this into four subcases: 

1. Let dxa{x,y,p) = dya{y,x,p) = 0, i.e., a = a{p). Then (1B.23P implies g{x,p) = x + g{p), 
which again yields a subcase of (jB.lSp . 

2. Le t dsa {x,y,p) = d^g{x,p) = 0, i.e., a{x,y,p) = a{y,p), g{x,p) = gi{p)x + go{p). 
Equation ()B.23P then reduces to (giip) — l)(y9ja(y, p) + a{y,p)) = 0. This first factor again 
yields a subcase of (jB.lSh but the second factor yields a new solution 

a = a{p)z/y, 7 = a{p){gi{p)x + c/2(p))/y, 9i{p) 7^ 1- (B.28) 

3. Next assume dxa{x,y,p) / 0, d'^g{x,p) = 0. Then we must in fact have dxg{x,p) = 1, 
leading back to (jB.lSh . 

4. Finally assume d^a{x,y,p) / 0,d'^g{x,p) ^ 0. Then dividing ()B.27p by their product and 
computing a x derivative yields the equation 

2(9.5 -1)^^5 = 3(925)2. (B.29) 

This has the solution 

g{x,p) = X + gi{p) + g2{p)/{x + gsip)). (B.30) 

Here we may assume 52 7^ 0, else we are back to a subcase of (lB.17p . Returning to ()B.23p we 
get from the coefficient of x two equations, 

ydya{x, y, p) + a{x, y, p) = 0, {x + 53)9ja(x, y, p) - a{x, y, p) = 0, 

which are solved by 

a{x,y,p) = a{p){x + g^{p))/y. 

This leads to the solution 
a{x,y,z,p) = a{p){x + gs{p))z/y, j{x,y,z,p) = a{p)[g2{p) + {x + g3{p)){x + gi{p))]/y. (B.31) 
Thus we have obtained four solutions, a = 0, (IB.lSp . f lB.28p . and (]B.3ip . 
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